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A new (algebraic) approximation scheme to find global solutions of two point boundary value 
problems of ordinary differential equations (ODE's) is presented. The method is applicable for 
both linear and nonlinear (coupled) ODE's whose solutions are analytic near one of the boundary 
points. It is based on replacing the original ODE's by a sequence of auxiliary first order polynomial 
ODE's with constant coefficients. The coefficients in the auxiliary ODE's are uniquely determined 
. from the local behaviour of the solution in the neighbourhood of one of the boundary points. To 

■ obtain the parameters of the global (connecting) solutions analytic at one of the boundary points, 

reduces to find the appropriate zeros of algebraic equations. The power of the method is illustrated 
by computing the approximate values of the "connecting parameters" for a number of nonlinear 
, ODE's arising in various problems in field theory. We treat in particular the static and rotationally 

Q ■ symmetric global vortex, the skyrmion, the Nielsen-Olesen vortex, as well as the 't Hooft-Polyakov 

V-^ ' magnetic monopole. The total energy of the skyrmion and of the monopole is also computed by 

the new method. We also consider some ODE's coming from the exact renormalization group. 
The ground state energy level of the anharmonic oscillator is also computed for arbitrary coupling 
CNj ' strengths with good precision. 

T— I . PACS numbers: 02.30Hq,02.30Mv,02.60Lj,11.27+d 
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It occurs quite often in Physics (but of course also in other areas of Science) that we have to solve (singular) 
two point boundary value problems (2p. BVP) associated to a system of ODE's. An important class of such BVP's 
for linear ODE's arises from eigenvalue problems of the stationary Schrodinger equation either in one dimension or 
reduced to an ODE with some (e.g. rotational) symmetry. Another large class of 2p. BVP's for nonlinear ODE's 
stems from the equations of motion of classical field theories reduced to ODE's (e.g. with some symmetries) and 
one could easily continue the list. We start with the example of the static, rotationally symmetric global vortex in 
a Ginzburg-Landau effective theory which has numerous applications ranging from condensed matter to cosmic 
strings, see 2, 3]. The field equation determining the vortex profile can be written as 

Qh. 

<D . In 2 

S"{r) + ~/'(r) + (1 - — )f(r) - f(r) = , (1) 
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where f(r) is a real function, /' = df / dr, and n £ Z corresponds to the vorticity. The physically interesting, globally 
regular solutions of Eq. ([T|) satisfy the following boundary conditions (BC): 

f(r^0) = k n r n + O(r n + 2 ), f(r - oo) = /«, = 1 , (2) 

and then a major problem of the 2p. BVP amounts to find the value(s) of the free parameter, k n , to ensure the BC 
of f(r) at r = oo. 

It is not too difficult to integrate numerically Eq. |T]) e.g. by the "shooting" method from r = to some large value of 
r and determine k n to some precision but it is considerably more difficult to obtain analytical results. The aim of this 
Letter is to present a new analytic procedure to approximate the value of k n for the connecting trajectory involving 
only algebraic steps. Its basic input is the power series expansion of the solution around a point of analyticity (typically 
at r = 0) and the BC at r — oo. With this input our method reduces the connection problem for k n to find the 
corresponding root of a polynomial equation. The method is conceptually very simple, it is easy to apply and moreover 
it yields good approximations for various ODE's. Our method is heuristic, we cannot put it on a mathematically 
rigorous footing as yet, nor can we precisely define the class of ODE's to which it is applicable. Nevertheless, it 
seems to us that with regard to its simplicity and its large applicability, the new method is of considerable interest for 
many applications (it yields with little effort good results for the connection parameters of many nonlinear ODE's, 
the energy levels of the quartic anharmonic oscillator for arbitrary values of the coupling, etc.). 

We illustrate our method in detail on the example of the global vortex fl}. The first step is to introduce the 
following sequence of auxiliary first order polynomial (implicit) ODE's of the form 

F N (f, f) := f' N + G 1 (f)f N - 1 + ... + G N ^(f)f + G N (f) = 0, N =1,2... , (3) 
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where Gi(f) = J^j Gij ft is a polynomial in / with constant coefficients. Such implicit ODE's are not easy to handle 
in general, however, as it will be shown here, one can squeeze out some important information from the sequence 
{F N (/' , /) = 0} without having to solve them. The next step in our method is the determination of the unknown 
coefficients Gy in F N (/' , /). In order to do this we enforce that Eq. ((3]) be satisfied to the highest possible order in 
r using the power series expansion of the solution of Eq. ([1]) f(r), around r = 0. Having found the constants Gij this 
way, we can now impose the BC at r = oo for f(r) in Eq. (|3|) which amounts to 

J2G Nj (uy = o. (4) 
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Since the coefficients Gij depend on k n (they turn out to be rational functions), Eq. (QJ represents a polynomial 
equation for k n . There is no a priori condition on the degree of the polynomials £?$(/), we have chosen to impose 
deg(Gi) < 2i. We remark that this restriction on the degree is known to be a necessary (but not sufficient) condition 
for the absence of movable branch points in Eq. |3]). Our main observation is that in the set of real roots, Sjv, of ||3J| one 
can find a root, rjv, which seems to converge to the value of the connection parameter and the corresponding trajectory, 
/Ar(r), of ||3J) yields a global approximation to the solution of the 2p. BVP of Eq. (JTJ). The main problem in our method 
is to identify the "good" root in Sn- At this point we also note that our method performs a kind of resummation 
from the local power series expansion through the auxiliary ODE's, nevertheless we see no obvious relation to more 
standard resummation techniques such as the Borel technique or Pade approximants, see the monograph Q for a nice 
review of these and other approximation techniques. 

We now show in detail how our method works for the simplest case N = 1, and for vorticity n = 1. Without loosing 
generality one can assume ki > 0. In order to obtain a nontrivial result for N = 1 we have to take in Eq. Q for 
Gi(f) a polynomial of degree two, i.e. the simplest auxiliary ODE is a Ricatti eq.: 

/' + G 10 + G n / + G 12 / 2 = 0. (5) 

Using the power series expansion of f(r) (by solving Eq. (JTJ)) f(r) — k%r — fcir 3 /8 + . . ., it is easy to obtain the 
coefficients Gu in Eq. ([5]): Gio = — fci, Gn = 0, G12 = 3/(8fci). Therefore the solution of Eq. (|4]) corresponding to 
the BC © is fci = a/3/8 = 0.612 . . . , which constitutes a reasonable first approximation for ki, (see Table I for the 
numerical value, fci nU m) and by solving Eq. (J5J) one obtains quite a good approximation for the vortex profile function, 

To improve upon this approximation one could keep N = 1 fixed and increase only the degree of Gi(/), however, 
since this scheme is not sufficiently general, (it works quite well in some, but not in all cases), we rather consider the 
next member, N = 2, in our auxiliary ODE sequence with deg(G^) = i. In this case we have to expand f(r) in Eq. (JTJ) 
up to order 5, and repeating the same procedure as for N — 1, one finds for the coefficients Gij: Gio = ki(80kf + l)/Di, 
Gn = 0, G 20 = -2kf(20kf + 7)/Di, G 2 i = 0, G 22 = 81/(8L>i), with D 1 = 13 - 40fcj\ Therefore Eq. © reduces to 

320fcf + 112*^-81 = 0, (6) 

whose positive real root is k\ = ( — 7 + \/454)/40 = 0.598..., which compares rather satisfactorily with fcinum, 
considering the simplicity of the calculations. At this stage it is natural to ask how this approximation changes 
keeping N — 2 fixed but increasing the degree of Gj to 2j. Repeating the computations for this case Eq. ((4]) yields 
1 + 368kf - 2400fcf + 3840^ = 0, which has two real positive roots, 0.587 . . . and 0.531 . . .. Now 0.587 . . ., is a better 
approximation than the previous one, however there is also a second "spurious" root, 0.531 . . ., which needs to be 
excluded. The appearance of such "spurious" roots for increasing N and for increasing degrees of Gj is a general 
feature and it is a main drawback of the method. In practice, however, this does not necessarily causes too serious 
problems , since one can easily follow the "good" roots by continuity and checking their stability against changing N 
and the degrees of Gj. In the following, unless indicated otherwise, we shall choose the degree of the polynomials Gj 
to be j, this fixes the global degree of F N (f, f) to be N. In this case there are altogether (N + l)(N + 2)/2 constants 
Gij to be determined, and the order of the power series in r to be used has to be chosen accordingly. We present the 
approximate values of the connection parameter k n for n — 1, . . . , 4 up to N — 8 in Table I. Remarkably for n = 1 
the N — 8 approximation yields 8 correct digits compared with fcinum- For increasing values of n the order of the 
series in r must be also increased for a given degree of F N (/' , /) and consequently the degree of the polynomial in k n 
is greater. We stress that all computations are analytic, (they have been performed on a standard desktop computer 
using Mathematica 5.2), except to find the roots of the corresponding polynomial where numerical methods became 
inevitable. 

Next we show that our method also yields good results for an important eigenvalue problem, the determination of 
the energy levels of the quartic anharmonic oscillator in 1 dimension. Using dimensionless variables the Schrodinger 
equation can be written as: 



f"(x) + {2E-l3x 4 -x 2 )f(x) =0, 



(7) 
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TABLE I: Convergence of the approximants for the connection parameters of the global vortex for n = 1, 2, 3, 4 



N 


fcl 


k 2 


fc 3 


&4 


3 


0.585 








4 


0.5831 




0.021 




5 


0.58315 


0.1527 


0.025 




6 


0.583190 


0.1529 


0.0264 


0.0028 


7 


0.5831894 


0.15310 


0.026183 


0.0034 


8 


0.5831894936 


0.15309 


0.026185 


0.0033 


h 


0.5831894959 


0.1530991029 


0.02618342072 


0.00332717340 



where (3 is the coupling parameter and E is the energy eigenvalue [13| . The 2p. BPV for the ground state wave- 
function (which is even) can be put in the form f(x = 0) = 1, f(x — ► oo) = 0. For N = 3 one obtains the following 
algebraic eq. for the eigenvalue E 

- 6408£ 6 + 12960/3S 5 + 2356£ 4 - 2976/3£ 3 - 2(1440/3 2 + 133)£ 2 + 168(3E + 25 = . (8) 

The harmonic oscillator corresponds to (3 = 0, whose exact ground state energy is E = 1/2. For (3 — there is a 
single positive real root of Eq. ([8]), 0.5166 . . ., and the choice of the roots for N — 3, E3, in Table HT1 has been done by 
following this root as (3 has been varied. We summarize our results for N = 3 , 4 , 5 and 10 in Table [Til where we also 
compare them with the known ones in the literature § denoted as E v . One can thus see that our method yields very 



p 


E 3 


E 4 


E 5 


-Bio 


E v 


0.2 


0.51 


0.562 


0.5598 


0.5591455 


0.5591463 


1 


0.78 


0.704 


0.6984 


0.6961795 


0.6961758 


2 


0.90 


0.813 


0.8065 


0.8037773 


0.8037707 


4 


1.08 


0.965 


0.9548 


0.9515767 


0.9515685 


100 


2.85 


2.56 


2.502 


2.4983125 


2.4997088 


400 


4.49 


4.02 


3.931 


3.930989 


3.9309313 


2000 


7.65 


6.86 


6.692 


6.694321 


6.6942209 


40000 


20.7 


18.6 


18.13 


18.13751 


18.137229 


2 x 10 B 


76.3 


68.4 


66.76 






2 x 10 9 


763 


684 


667.6 







TABLE II: Convergence for the ground state energy of the anharmonic oscillator, En for N = 3 , 4 , 5 and 10 

good approximate values for the ground state energy of the anharmonic oscillator for all values of the coupling. We 
remark that some of the other real roots of the polynomial in E are related to the energy levels of the excited states 
and by our method one can also obtain approximate values for them, but we will not elaborate on this point here. 

Let us next present here some results on the simplest "skyrmion" solution, which is of considerable interest as a good 
approximation for baryons, we refer to the recent monograph [3j for details. The pertinent ODE for the spherically 
symmetric skyrmion field can be written as: 

(r 2 /')' + 2/" sin 2 / + sin(2/) [/' 2 - 1 - (sin 2 / )/r 2 ] = , (9) 

together with the BC's f(r) = n + kr + C(r 3 ), / — > for r — > 00. In this case when applying our method we have 
found that the choice deg(Gi(/)) = 2i yields significantly better results than deg(Gi(/)) = i. This way we find for 
the connection parameter: k — —2.084 for N = 2; k = —1.996 for N = 3 and k = —2.003 for N — 4, whereas 
&num = —2.007. The agreement is quite satisfactory taking into account the relatively low degree of the auxiliary 
equations ©. 

An important physical quantity is the total energy of such localized solutions. For example the energy of the 
skyrmion in a ball of radius R is given as E(R) — J Q £dr, where £ is the energy density: 

£ = [r 2 f 2 + 2(1 + f 2 ) sin 2 / + (sin 4 /)/r 2 ]/(37r) , (10) 

and the total energy is then E(oo). The direct way to compute approximate values for the energy, i.e. finding the 
corresponding solutions of the auxiliary ODE ([3]) first and then evaluating E(R) for them would be quite difficult 
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without resorting to numerics. We can apply a slight variation of our method, however, to obtain approximate values 
for the total energy once the value of the connection parameter is known. To do this we consider the sequence 
{F N (£,E) — 0}. Since the power series expansion of the energy, E(r), is completely fixed for any given value of 
the connection parameter, k, all the coefficients, Gy, will be determined. Then we solve the corresponding algebraic 
equation ((4]) for the unknown value of E(oo). This way for N = 9 we have obtained E(oo) — 1.22, which seems to be 
quite good when compared to the numerical value, E num = 1.23. 

We consider now some ODE's originating from a completely different problem, the fixed point equation of Wilson's 
exact renormalisation group (RG). The RG equation for scalar field theories in the local potential approximation can 
be written as Q: 

2f"(x) - if(x)f'(x) - 5xf(x) + f(x) = 0, (11) 

where f(x) = V'{x) — x, and V(x) is the potential. The pertinent solution of (llip is an odd function of x and for 
x — > f(x) — kx + 0(x 3 ). For large values of x f(x — > oo) — > ax 1 / 5 where a is a constant. Since in the present 
case foo — oo, /" — » and /' — > for x — > oo, it is rather natural to slightly modify the method by considering 
the auxiliary ODE's © for (/",/'), i.e. F N (f",f) = 0. Then proceeding exactly as before we obtain a polynomial 
equation in k, and we find with a good convergence k = —1.22859876 for N = 6. This agrees quite well with the 
numerical value fc num = -1.22859820 Q. 

Let us consider yet another example, the Wegner-Houghton's fixed point equation in the local potential approxi- 
mation m, 

2 ln(l + v"(x)) + 6v(x) - xv'(x) = . (12) 

The change of dependent variable, f(x) = v'(x), gives the simpler differential equation 

2f"(x) + [1 + f(x)} [5f(x) - xf(x)} = . (13) 

The solution of interest satisfies the following BC's: f(x — > 0) = kx + 0(x 3 ) and fix — > oo) — > ax 5 where a is a 
constant. The connection parameter, k, obtained numerically is fc num = —0.461533 [6]. The problem is similar to 
the precedent one but we use now higher derivatives, F N (f^ 7 ',f^') = 0. Proceeding in the same way that in the 
previous case we obtain k = —0.46144... for N = 6 which agrees with the numerical value, although not so well as in 
the previous case, probably due to the use of derivatives of higher order. 

We now show that our method can be generalized in a very simple way for a system of M ODE's for the set of 
unknowns {/ m (^)}i m = 1, ■ • ■ , M. We introduce for each unknown function a first order auxiliary implicit ODE 
{Fm m {f'rm fm) = 0}, where F^ m {f' m , f m ) is a polynomial of degree N m in f m (c.f. Eq. ([3|)). The constant coefficients 
in F^™ (fm> fm) are determined by demanding that {Fm m ifm>fm) = 0} be satisfied to the highest possible order 
in the power series solutions of f m at the origin, say. Then proceedings exactly as for the case of a single unknown 
we impose the BC at infinity and obtain a system of algebraic equations of the form J^j ^N m j(fm(°°)) = f° r the 
connection parameters. As a concrete illustration we shall consider the field eqs. of the static, rotationally symmetric, 
gauged vortex of Nielsen- Olesen Q and those of the 't Hooft-Polyakov magnetic monopole (Toj . 

The differential equations determining the cylindrically symmetric magnetic potential resp. scalar fields, a, / of the 
Nielsen-Olesen vortex with a single magnetic flux quantum can be written as: 

r(rf)' - / [(1 - af - r 2 (3(\ - f 2 )] = , (14a) 
ra" -a' + 2r(l - a)/ 2 = , (14b) 

where j3 corresponds to the self-coupling of the scalar field. The BC's necessary to ensure regularity and finite energy 
at r = and at r = oo are: f = fi = c\r + 0(r 3 ), a = /2 = c^r 2 + 0(r 4 ), /(oo) = 1, a(oo) — 1. To implement our 
method we have chosen F^ 1 (f[, fi) resp. F^if^, h) to be polynomials of degree N\ resp. N2 — N\ + 1 in order to 
have roughly the same number of terms in the algebraic eqs. for ci, C2. Next we recall the differential equations for 
the static spherically symmetric magnetic monopoles in a spontaneously broken SU(2) gauge theory: 

(r 2 tt)' /i[2/f + y/3 2 (A 2 - 1)] - 0, (15a) 
r 2 / 2 '-/2[(/|-l)+r 2 / 2 ]-0, (15b) 

where fi and correspond to the Higgs and the gauge field, respectively. The regular BC's at r — resp. r = 00 
are /1 = d x r + C(r 3 ), f 2 = 1 - d 2 r 2 + 0(r 4 ), /i(oo) = 1 and / 2 (oo) = 0. We present in Table EI] our results 
for the connection parameters of both the vortex and of the magnetic monopole for N\ — 2> and N 2 — 4, together 
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Cl 


C2 


Clnum 


C2num 


di 


d 2 


G^lnum 


d'2num 





- 


- 


- 


- 


0.3329 


0.1523 


1/3 


1/6 


1 


0.8542 


0.5007 


0.8532 


0.5000 


0.7343 


0.3422 


0.7318 


0.3409 


2 


1.0996 


0.6169 


1.0993 


0.6166 


1.0692 


0.4501 


1.0683 


0.4491 


3 


1.2846 


0.6975 


1.2843 


0.6969 


1.4025 


0.5350 


1.4003 


0.5321 


4 


1.4393 


0.7609 


1.4387 


0.7597 


1.7406 


0.6100 


1.7339 


0.5997 


5 


1.5748 


0.8137 


1.5741 


0.8119 


2.0872 


0.6901 


2.0701 


0.6567 


6 


1.6972 


0.8595 


1.6962 


0.8569 


2.4421 


0.7834 


2.4091 


0.7060 


7 


1.8097 


0.9001 


1.8082 


0.8967 


2.8035 


0.8938 


2.7503 


0.7492 


10 


2.1052 


1.0011 


2.1024 


0.9944 


3.6269 


0.9678 


3.7850 


0.8542 



TABLE III: Connection parameters for the vortex and the magnetic monopole with N\ — 3, N2 — 4 



with the corresponding numerical values from Refs. [ill . Il2j . The agreement is good, taking into account that the 
computations involving two polynomials are more cumbersome, moreover for both for small and large values of /?, one 
needs to increase the values of N\ and N2 more and more. We recall here that for (3 = eqs. fll~5l) are analytically 
soluble, yielding d\ = 1/3, di = 1/6. We have also computed the total energy for (3 — 1 applying the same procedure 
as for the skyrmion, and we have obtained E(oo) = 1.2136 for N = 10 with deg(Gj) = 2j (c.f. E num = 1.237 [Hj]). 

In conclusion we have presented a new algebraic scheme to obtain the connection parameters, the energy eigenvalues 
and the total energy for a number of physically interesting two point boundary value problems associated with ODE's. 
The method is based on a sequence of auxiliary first order implicit polynomial ODE'S, which is determined from the 
series expansion of the solution to a given degree. Imposing the BC's for the solution of the auxiliary ODE's yields 
algebraic equations for the connection parameters. The computation of power series expansions and finding the roots 
of the algebraic equations are easily implemcntablc on symbolic formula manipulation systems. It seems to be an 
interesting problem to clarify the mathematical basis of our method. 

Acknowledgements We would like thank Claude Bervillier for pointing out that the fixed point equations of the 
RG serve as interesting illustrations of our method and for his interest. We also thank him and Peter Breitenlohner 
for providing us various numerical results. One of us (B.B.) had also some interesting discussions with Stam Nicolis. 



[1] V.L. Ginzburg and L.P. Pitaevskii, Sov. Phys. JETP 34, 858 (1958); L. P. Pitaevskii, Sov. Phys. JETP 13, 451 (1961). 
[2] A. Vilenkin and E.P.S. Shellard, Cosmic Strings and Other Topological Defects, Cambridge University Press, (1994). 
[3] N. Manton and P. Sutcliffe, Topological Solitons, Cambridge University Press, (2004). 

[4] CM. Bender and S.A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, McGraw Hill, (1978). 
[5] R. Karrlein and H. Kleinert Phys. Lett. A 187, 133, (1994). 
[6] C. Bagnuls and C. Bervillier, Phys Rep. 348, 91 (2001). 
[7] C. Bervillier, private communication. 

[8] A. Hasenfratz and P. Hasenfratz, Nucl. Phys. B 270 [FS16], 687 (1986). 

[9] H.B. Nielsen and P. Olesen, Nucl. Phys. B 61, 45 (1973). 
[10] G. »t Hooft, Nucl. Phys. B 79, 276 (1974); A. M. Polyakov, JETP Lett. 20, 194 (1974); 
[11] P. Forgacs, S. Reuillon and M.S. Volkov Nucl. Phys. B 751, 390 (2006). 
[12] P. Forgacs, N. Obadia, and S. Reuillon, Phys. Rev. D 71, 035002-1 (2005). 
[13] Note that our j3 corresponds to g/2 of Ref. [f| 



